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f~>) , The existence of the spectral break around ~ 3 x 10 eV in the cosmic ray spectrum 

^v^j . (referred to as the 'knee') is one of the biggest questions in cosmic ray astrophysics. At the 

same time, the origin of cosmic rays above the knee energies (between 10 15 and 10 18 eV) is 
' also still unsettled. In this paper, we investigate how the hypothetical extragalactic CRs after 

r ^ , modulated by the galactic wind contribute to the knee in the CR spectrum. We numerically 

calculate the modulated energy spectrum of the hypothetical cosmic rays coming into the 
' galaxy from just outside of the "galactic sphere" where the galactic wind terminates. We 

show that the observed knee structure is reproduced well by a superposition of the modulated 
, component and the galactic cosmic rays originating in supernova remnants. 

>■ 

<N ' 

^r) [ §1. Introduction 

Observations have revealed cosmic rays (CRs) spread over 11 decades of energy 
from ~ 10 9 eV to ~ 3 x 10 20 eV. The spectrum is often fitted by power laws of 
energy with the index of ~ —2.7 below ~ 3x 10 15 eV and ~ —3.0 above, respectively. 
This break in the spectrum is referred to as the 'knee' in the spectrum. The knee 
has been recognized as an important structure of the energy spectrum, which may 
provide constraints on the acceleration and propagation of CRs. 

CRs with energies below the knee have been believed to be originated in su- 
pernova remnants (SNRs) in our Galaxy from general arguments about energetics 
and the diffusive particle acceleration mechanism in shocks. 1 **' 2 ) This argument 
of the CR origin is corroborated by recent observations of X-rays (by ASCA) and 
TeV gamma-rays (by CANGAROO and HESS) from two SNRs, SN1006 3 )' 4 ) and 
RX J1713.7— 3946. 5 - H9 - ) These observations have revealed that electrons and protons 
are accelerated to energies of ~ 100 TeV in the SNRs. *) However this maximum 
energy does not reach the knee energy of ~ 3 x 10 15 eV. A hypothesis blessed by 
long tradition explains the break in the spectrum at the knee by more rapid escape 
from the Galaxy of more energetic particles above the knee. 11 ) Alternative interpre- 
tations of the knee structure are to invoke one or a few extra components of galactic 
origin which dominate above ~ 1 PeV. 12 ^ 15 ^ Yet the origin of CRs above the knee 
is still not settled. 

Recently the existence of diffuse CR electrons in the intergalactic space has been 
suggested from the results of extreme-ultraviolet and high energy X-ray observation 
of clusters of galaxies 16 -* and from the observation of the diffuse cosmic gamma-ray 
background. 17 - )> 18 ) If nuclear components with energies extended well above ~ 1 PeV 
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*^ For SN1006, there is discrepancy between the results of CANGAROO-I and HESS; an upper 
limit reported by HESS is below the flux detected by CANGAROO-I. 10) 
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also exist together with the diffuse electrons around our Galaxy, these components 
modulated by the galactic wind might be directly observable at the earth. In this 
paper, we numerically examine such a possibility, and discuss their implications for 
the origin of the knee. 



§2. Numerical Simulations 



We postulate the existence of hypothetical CRs just outside of the "galactic 
sphere" where the galactic wind 19 ) terminates. The energy spectrum of these CRs 
is assumed to be the same as the spectrum of the CRs observed at the earth with 
energies higher than the knee region but extrapolated to much lower energy range; 
namely the spectrum is proportional to E~ 3 , where E is the total energy of a particle. 

These CRs may diffuse into inner region of the galactic sphere against the ex- 
panding galactic wind. We examine how the spectrum of these CRs should be 
modulated during this propagation process. The transport of CRs is described by 
the Fokker-Planck equation for the spherical symetric case 20 ) -22 ) 

dt r 2 dr dr dr r 2 dr 3 dp ' 

where / is the phase space distribution function, t is the time, r is the radial distance, 
V is the speed of galactic wind, p is the particle momentum, and k is the diffusion 
coefficient for radial propagation. Here we neglect energy change processes other 
than adiabatic losses. 

It is known that Eq.Q is equivalent to the coupled stochastic differential equa- 
tions (SDEs). The SDEs equivalent to Eq.Q are written using new quantities 
u = \n(p/mc) (where m is the particle mass and c is the speed of light) as 

dr = (V + — )dt + V2ndW r , (2) 
r 

and 

2V 

du = dt, (3) 

3r 

where dW r is a Wiener process given by the Gaussian distribution, P(dW r ) = 
(2-Kdt)~ l l 2 exp(— dW 2 /2dt). Here we assume that V does not depend on r. If we 
integrate the SDEs "backward in time", 23 )' 24 ) we can obtain a probability function 
F(p, R\po, r ) which is necessary to calculate the modulated energy spectrum at an 
arbitrary point ro from the Galactic Center. The probability function describes the 
probability of which a particle observed with momentum po at ro would have had 
momentum p at the boundary of the galactic sphere R. Once we have calculated 
F(p, R\po, ro), the modulated energy spectrum f ro (po) at a point ro is calculated 
with the energy spectrum fn(p) at R as, 

fr (Po) = J f R (p)F(p,R\p ,r )dp. (4) 

This method has been applied successfully to the investigation of the solar mod- 
ulation phenomena of the galactic CRs in the heliosphere. 23 - 1 ' 24 ) The momentum 
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spectrum of the hypothetical CRs at the boundary R is Jr oc p~ 5 , if we assume the 
spectrum is proportional to as mentioned earlier. 

We integrated numerically Eqs.© and © assuming the diffusion coefficient n 

as 



where r] is the ratio of the mean free path of the particle to the Larmor radius and k-q 
is the Bohm diffusion coefficient, k-q = Ec/(3ZeB) (where E, Z, and B are the total 
energy of particle, the atomic number and the magnetic field intensity, respectively). 

The resultant modulated spectrum at an arbitrary point ro depends on five pa- 
rameters, 77, Z, R, V, and B, where these parameters are unknown to us. According 
to the study of the galactic MHD wind driven by CRs in a rotating galaxy, 25 ) it is 
expected that our Galaxy is surrounded by a large galactic wind halo with a scale of 
the order of R ~ 100 kpc. Also, the resultant velocity profile reveals that it increases 
almost linearly with distance from the galactic disk up to 20 kpc and then becomes 
a constant with about ~ 300 km s -1 up to ~ 100 kpc. When the extragalactic 
CRs propagate into inner region of the galactic sphere, it will spend long times in 
the outer region of the galactic sphere rather than the inner region like around the 
galactic disk. Therefore we should select R ~ 100 kpc and V ~ 300 km s -1 in our 
calculation. 



3.1. Galactic Modulated Spectrum 

Fig.^a) shows the calculated differential energy spectra of protons (Z=l) at the 
earth (ro = 8.5kpc) as a function of total energy E. The solid, dashed, and dotted 
lines are the differential intensity with r/=1000, 100, and 10, respectively. Here we 
assume V = 300 km sec -1 , R = 100 kpc, and B = 1/zG. The straight line indicates 
the assumed unmodulated spectrum at the boundary of the galactic sphere. The 
break point of the spectrum should be compared with the knee. We can find that 
the spectrum is break around the knee energy. We also find that if r\ increases by 
some factor, the break point is shifted to the lower energy by the same factor. We 
also calculated the modulated energy spectra of hypothetical nuclear components 
instead of protons. The results are shown in Fig. ^b), where each lines are the 
differential intensity of components labeled in the figure at the earth with rj = 1000 
as an example. We find that the break point of the nuclear components is shifted 
to the higher energy by a factor of Z compared with that for protons. It may also 
be worth remarking how is the break point shifted when the other parameters, B, 
R, and V are changed. The resultant spectra are shown in Fig. ^c)-(e). Here if we 
define the breaking energy inbreak as the energy in which the modulated spectrum 
shown in Fig. ^ becomes a maximum, it is approximately shown as 




§3. Results 



E hieSLk {Z, V ,B,R,V) ~ 6 x 10 15 Zrj^BAV* eV, 



(6) 
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Fig. 1. Expected energy spectra of the hypothetical CRs modulated by the galactic wind. The 
straight line indicates the differential energy spectrum of the hypothetical CRs with Z at the 
boundary of the galactic sphere 7? which is a power law in total energy with a spectral index of 
—3.0. (a)the differential intensity of protons (Z=l) at the earth (ro = 8.5 kpc) for R — 100 kpc, 
V = 300 km sec -1 , and B = 1/xG with r] = 1000(solid line), 100(dashed line), and 10(dotted 
line), respectively. (b)Same as (a) but with five values of Z for r\ — 1000, R = 100 kpc, V = 300 
km sec -1 , and B = l^iG. (c)Same as (a) but with four values of R for r\ = 100, Z=l, R = 100 
kpc, V = 300 km sec -1 , and B — lfiG. (d)Same as (a) but with four values of V for 77 = 100, 
Z=l, R = 100 kpc, and B = 1/iG. (e)Same as (a) but with four values of B for r\ — 100, Z =1, 
R = 100 kpc, and V = 300 km sec -1 , (f) Same as (a) but at four values of an arbitrary point 
r for r? = 100, Z=l, R = 100 kpc, V = 300 km sec -1 , and B = 1/iG. 
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Kinetic Energy T (GeV/nucleon) 

Fig. 2. Differential energy spectra obtained by various direct measurements as a function of kinetic 
energy per nucleon and the best fit curves for each element. The data were adopted from 
Refs. 26)-31). 



where ^=r//(100), B*=B /(1/zG), i?*=i?/(100kpc), and K=V/(300 km s" 1 ), respec- 
tively. This results indicate clearly that we can reproduce the knee by choosing an 
appropriate values for these parameters. 

In addition, in order to investigate how the spectrum is modulated in arbitrary 
points other than the earth, we calculated the modulated spectra at four values of 
tq. The results are shown in Fig. ^f). We found that if ro approaches the boundary 
R, the break point of the spectrum shifts to low energy more and the modulated 
spectrum diminishes more rapidly as energies go down at the same time. This result 
indicates that the hypothetical extragalactic CRs with energy lower than at least 
TeV cannot penetrate into the inner region of our galactic sphere due to the effect 
of modulation by the galactic wind. 

3.2. Model Spectrum (Two Components Model) 

In this section, to demonstrate how well our model reproduces the observed 
all-particle spectrum near the knee region, we consider a simple model which is 
mentioned below. We assume that the all-particle spectrum observed at the earth 
-^totai(-S') is superposition of the two components, namely the modulated extragalac- 
tic component F mo( ^ u \(E) and the component originated in SNRs in our Galaxy 
Fsnr(E) as 



-^total — FSNR(E) + -^modul 



(7) 
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and Fsnr(-E') and -FmoduiC-^) are represented as 

^snr(^) = ^/ s Z nr(^) (8) 

z 

Fmodul(E) = ^2 f^ odu i(E), (9) 

z 

where f§^(E) and f^ ^ u \(E) are each nuclear components, respectively. Now we 
consider that Fs^r(E) consists of the sum of components after propagated in in- 
tersteller space from each of SNRs in our Galaxy, in which particles are accelerated 
to a power law by shock acceleration. Thus we simply assume that each nuclear 
component observed at the earth /§nr(E) is represented by a power law with each 
spectral index of a as 

fi NR (E) oc E- a exp(-E/(ZE m&x )), (10) 

where -E max is the maximum energy of protons which may come from the maximum 
energy attained by protons accelerated in SNRs. Here we reexamined published data 
of CRs which are believed to originate in SNRs in our Galaxy as shown in Fig. [2] 
as a function of kinetic energy per nucleon T. Here we fit with f§^j>(T) oc T~ a exp 
(— T/(ZE max / A)) because the published data is represented as a function of T. We 
assumed in this fitting E max = 500 TeV for example, and also assumed that (Z, A) 
are (7, 14) for CNO-group and (12, 24) for NeMgSi-group, respectively. We define 
the sum of these components shown in Fig. [2] as SNR component Fsnr(E). 

Fig. El show the resultant model spectrum F to t[d(E) around the knee for two 
example cases with the all-particle data obtained by various experiments. Fig. E^a) 
shows a model spectrum when we adopt F modu \(E) for the modulated proton spec- 
trum /modu^-E 1 ) calculated with rj = 250, R = 100 kpc, V = 300 km sec -1 , and 
B = 1/xG. As clearly seen in Fig. Efa), we find our model reproduces the observed 
spectrum around the knee fairly well. We emphasize that we can replace the mod- 
ulated proton spectrum / mo( j u i(-^) with another nuclear component f^ odn \{E) cal- 
culated with the values according to the roles shown in Eq.@. For example, if we 
adopt rj = 500, R = 100 kpc, V = 300 km sec -1 , and B = 0.077/uG, we can replace 
/modul ( E ) with the modulated iron spectrum f^ oAul (E). 

On the other hand, we don't know the chemical composition of the hypothetical 
extragalactic CRs. As an extreme example, we assume compound F mo d n i(E) which 
consists of the same spectrum of / modu] (£), & ul (E), f™° ul (E), Gf(E), and 
/moduli) a * ^ ne boundary. The result is shown in Fig.|3fb), where we adopt r\ = 250, 
R = 100 kpc, V = 300 km sec -1 , and B = 1/iG, and each abundance ratio at 
the boundary is five times smaller than that for the case in Fig. Efa). We can 
also reproduce a spectrum by using the model with compound of F mo d u \(E). We 
demonstrated our model for a simplistic case of spherical symmetric geometry with 
specific forms of K. 

3.3. Mean Mass of CRs around the Knee 

On inspection of Fig. |21 we also find the mean mass of the CRs around the 
knee should change with energy in complicated manner, because SNR component 
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Fig. 3. Observed all-particle spectrum and the model spectrum around the knee versus total energy 
of particle E. The observed data were adopted from Refs. 27),30)-34). The thick solid line 
indicates the model all-particle spectrum Ftotn\{E) which is the superposition of two components 
shown by two solid lines; SNR component Fsnr(-E) and modulated extragalactic component 
F mo d u i(E). The SNR component is the sum of all nuclear components presented in Fig. |2]but 
presented here as functions of different energy scale, (a) Single component model for F mo dui{E). 
Modulated component is calculated for proton /,^ odul (£) with r\ — 250 and with the same 
values for R, V, and B as in Fig.Q^a). Straight dotted line indicates the energy spectrum of the 
hypothetical CRs at the boundary R as ~ 4.3 x 10 20 E~ 3 cm -2 sec -1 str -1 eV -1 . Notice that 
we can replace f^ odu i(E) with another nuclear component fmadvd(E) (see in the text in details), 
(b) Composite model for -Fmodui(-E). Modulated component consists of the sum of five chemical 
comonents (f» odnl (E), f^ odul (E), /£od°ul(£), /£3T0B), and f^ odul (E)). Each components is 
calculated with r\ = 250, R = 100 kpc, V = 300 km sec , and B = 1/iG, where each spectrum 
assumed at the boundary R is the same (shown as straight dotted line) and is five times smaller 
than that for the case in Fig. |^a). In both figures, the numerical factor in n was chosen so that 
the model spectrum reproduces best the observed spectrum. 



diminishes gradually depending on its charge and the modulated component creeps 
in as energies go up. We have to estimate the expected mean mass to be compared 
with measurements by assuming the chemical composition of the hypothetical CRs, 
because we do not know it. By using the results of our model spectrum shown in 
Fig. |3J we calculated mean mass of the CRs. The results are shown in Fig. 21 as a 
function of total energy E with the mean mass determined by direct measurements 
in regions below the knee. Dashed and solid lines indicate the case for a model 
with F modul (E) = C odul (E) and with F modul (E) = f^ odnl {E) in Fig. Ha), and 
dotted line with a composite abundance for F modn \(E) in Fig. E{b), respectively. 
Experimental determinations of the mean mass in the region above the knee are, 
however, scattered in rather wide range for various cause depending on the methods 
in ground-based measurements and also on the interaction models of hadrons in the 
data analysis. 35 )~ 37 ) Our model predicts the mean mass above ~ 50 PeV would tend 
to that of the hypothetical extragalactic CRs as inferred from Figs. |3] and 01 Our 
model should be testified by future experiments in the energy range much higher 
than the knee. 
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§4. Discussion 

4.1. Energetics of the hypothetical CRs 

Heretofore we have simply assumed the existence of the hypothetical extragalac- 
tic CRs without any specification of their origin. We should address some comments 
on their origin including energetics. As an extreme case, if the hypothetical CRs 
pervade the intergalactic space uniformly and the spectrum is extended down to 
their rest mass energy as ~ 4.3 x 10 20 £- 3 cm 2 sec 1 str 1 eV 1 which we have 
assumed in Fig. Eta) (straight dotted line), their energy density is estimated to be 
as high as ~ 294 eV cm -3 when we assume they are protons. The density parameter 
corresponding to this energy density is J2/i 2 ~ 0.057 which is much larger than that 
of total baryon in the universe of QhjQ ~ 0.020 invoked by the Big Bang Nucle- 
osynthesis. 38 ) Therefore it is natural to consider that either the spectrum becomes 
harder and / or there is a cut off by other reasons in the energy region lower than the 
knee, and/or the hypothetical CRs are confined in local regions surrounding galaxies 
as well as our Galaxy and/or cluster of galaxies. Actually if we assume that the 
hypothetical CR spectrum becomes harder with the index of 2 in the energy range 
lower than 10 14,5 eV, then the expected energy density is estimated to be as small 
7.9 x 10~ 3 eV cm without any influence on the resultant modulated spectrum. 
Yet there might remain a possibility that some portion of the density parameter 
comes from these invisible hypothetical CRs. 

4.2. Hypothetical CRs in SMC 

Measurements with the EGRET telescope show that the density of CRs inside 
the Small Magellanic Cloud (SMC) is several times less than in our Galaxy. 39 ) This 
fact indicates simply that GeV CRs may originate in each galaxy, possibly in SNRs. 
This results seem to be incompatible with the existence of the hypothetical extra- 
galactic CRs which we introduced here. However if SMC have the galactic wind like 
our galaxy, the intensity of hypothetical CRs around GeV region inside SMC should 
be suppressed due to the modulation effect of SMC itself. Actually the galactic wind 
may have collided with that of SMC, LMC, etc., and it may have the complicated 
structure in which the shock waves are produced on the boundary. In this paper we 
do not address its detail. To investigate the behaviour of the hypothetical CRs in 
details, we need further simulation study for complex model using three-dimensional 
code including a few galactic sphere (for our Galaxy, LMC, SMC, etc.). 

4.3. Possible Origin of Hypothetical CRs 

There are some scenarios to support the existence of the hypothetical extra- 
galactic CRs which we introduced here. Volk & Atoyan 40 ) discussed the existence 
of non-thermal hadronic components in the cluster of galaxy in terms of early star- 
bursts and magnetic field generation in galaxy clusters. They suggest that relativistic 
particles would be confined in galaxy clusters over times longer than the age of the 
universe. Another interesting scenario is that it is possible to accelerate CR protons 
up to ~ 10 18 eV by cosmic shock waves which are generated by large scale structure 
formation 41 ^' 42 ) and by mergers between clusters of galaxies. 43 ) Miniati et al. 44 ) 
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Fig. 4. Mean mass of CRs around the knee as a function of total energy of particle. Symbols indi- 
cate the values obtained by direct measurements. 30 - 1, 31 - ) Dashed line; -Ft t a i(-B) = Jsnr(S) 



+ /moduli) ( shown in Fig- H a ))- Solid line; F total (E) = F SNR (E) + f^ odul (E) (shown 
in Fig. Ha)). Dotted line; F total (£) = F SNR (E) + fLam(E) + f"Z dal (E) + /SSSi(S) + 
fZaf L (E) + fZlUE) (shown in Fig.^b)). 



have carried out a computational study of production of CR protons at cosmologi- 
cal shocks and pointed out that a significant fraction of the total energy associated 
with baryons inside a cluster could be stored in CRs as a consequence of diffusive 
particle acceleration at structure formation shocks. In the future, GeV— TeV 7 ray 
observation of the violent cluster merger event shown in the evidence for shock accel- 
eration 45 )' 46 ) might directly lead to the existence of shock accelerated extragalactic 
CRs and might clear the relation between the extragalactic CRs and the CRs above 
the knee energy. 



We show that the galactic wind can provide a cutoff in about the right energy 
range lower than that around the knee. The all-particle spectrum of CRs around the 
knee is explained by the superposition of CRs originated in SNRs in our Galaxy and 
the hypothetical extragalactic CRs modulated by the galactic wind. The CRs above 
the knee region may be survivors of the extragalactic CRs in the battle against the 
outflowing galactic wind. 
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